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asymptotic  expansion  method  was  used  to  derive  a  relationship  between  the 
parameters  characterizing  the  structural  nonlinearity  and  the  amplitude 
and  frequency  of  the  limit  cycle  response. \  First  and  second  order 
perturbation  solutions  were  obtained  to  moael  the  "effective"  or 
linearized  system  parameters  governing  th^/nonl inear  response  of  the 
undamped  system. 


^-VThe  results  of  this  investigation  show  that  the  asymptotic  solutions 
accurately  predict  the  stationary  limit  cycle  behavior  when  compared  with 


numerical  simulation  and  describing-function  analyses  for  the  non- 
linearities  considered.  The  influence  of  higher  harmonics  on  the 


predicted  limit  cycle  response  were  also  observed  when  higher  order 
perturbation  solutions  were  obtained. 

The  aeroelastic  response  of  a  baseline  aerodynamic  surface  was 
investigated  using  the  asymptotic  expansion  approach.  Flutter  results 
were  obtained  for  the  effective  system  for  both  rigid  and  flexible 
representations  of  the  surface.  .Steady  state  aerodynamic  loading  was 
assumed  in  computation  of  the>Tutter  results.  These  results  were  then 
utilized,  in  conjunction  gitn  the  asymptotic  solutions,  to  investigate 
the  interrelationshijj^-between  the  magnitude  of  the  nonlinearity,  flight 
condition  andUtheTlmit  cycle  response.  Examples  and  results  of  the 
appli^atrOTTof  the  developed  solutions  are  presented. 

VThis  study  demonstrates  the  applicability  of  the  asymptotic 


expansion  method  in  accounting  for  the  influence  of  structural  non- 

linearities  in  the  limit  cycle  analysis  of  aerodynamic  surfaces. _ 

Understanding  the  feffects  of  structural, non linearities  on  the  dynamic 
response  of  systems  important  when  considerations  of  stabi 1 i ti/\ 
response  and  fatigue  life  influence  the  design  of  the  surface  sijbport 
structure  and/or  actuators.  \ 
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1.0  INTRODUCTION 


Defining  the  flutter  and  divergence  xcharacteri sties  of  aerodynamic  sur¬ 
faces  is  a  basic  requirement  to  assure  structural  and  performance  integrity  of 
a  given  design  in  its  operational  environment.  For  systems  which  have  struc¬ 
tural  nonlinearities  another  mode  of  aeroelastic  response,  limit  cycle  oscil¬ 
lation,  may  be  present.  A  comparison  between  limit  cycle  response  and  clas¬ 
sical  flutter  and  divergence  is  illustrated  in  Figure  1.  The  limit  cycle 
response.  Figure  1(c),  is  defined  as  a  constant  amplitude  steady  state  oscil¬ 
lation  whereas  the  divergence.  Figure  1(a),  and  flutter.  Figure  1(b),  are 
unstable  motions  with  increasing  amplitude.  The  importance  of  the  limit  cycle 
response  is  the  potential  of  these  oscillations  to  occur  within  the  flutter 
and  divergence  flight  envelope.  Frequently  real  hardware  designs  do  have  non- 
linearities  in  the  surface  support  structure  and/or  actuators  as  a  result  of 
manufacturing  tolerances  and/or  freeplay.  When  these  nonlinearities  exist, 
the  classical  assumption  of  a  linear  force-displacement  relationship  is  no 
longer  justified  and  an  understanding  of  the  nonlinear  effect  on  the  dynamic 
behavior  is  required  to  evaluate  the  system  response.  In  this  study  effort 
analysis  procedures  were  developed  to  characterize  the  limit  cycle  response  of 
aerosurfaces  with  discrete  structural  nonlinearities. 

The  effects  of  structural  nonlinearities  on  aerodynamic  surface  response 
has  been  studied  both  analytically  and  experimentally.  Reference  1  through  5. 
In  these  studies  several  nonlinearities  that  are  typically  encountered  in 
aerodynamic  surface  designs  were  considered.  In  the  analytical  studies  of 
References  1  through  4  the  method  of  harmonic  balance  or  describing-function 
approach  was  used  to  characterize  the  nonlinear  behavior  of  an  aerodynamic 
surface  with  root  freeplay  nonl inearities  of  the  type  shown  in  Figure  2. 

These  nonlinearities  are  representative  of  a  deadband  or  "slop"  in  the  root 
support  structure  with  and  without  a  linear  preload.  The  resulting  behavior 
is  such  that  the  stiffness  and  force  developed  in  the  adjacent  members  4s  a 
nonlinear  function  of  amplitude.  The  describing-function  approach  uses  a  one 
term  Fourier  Series  expansion  of  the  force  to  account  for  the  effect  of  this 
nonlinear  stiffness  on  aerosurface  response.  This  method  gave  satisfactory 
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(c)  Limit  Cycle  Oscillation 


FIGURE  1  MODES  OF  AEROSURFACE  AEROELASTIC  RESPONSE 


Dtietacanwnt 


(a)  FREEPLAY  NONLINEARITY 


(b)  PRELOAD  NONLINEARITY 


FIGURE  2  TYPES  OF  STRUCTURAL  NONLINEARITIES 


results  when  the  amplitude  of  motion  was  greater  than  the  magnitude  of  the 
freeplay,  Figure  2(a),  or  freeplay  plus  preload.  Figure  2(b).  However,  it  was 
pointed  out  References, 1  and  2,  that  when  the  amplitude  of  motion  is  approxi¬ 
mately  equal  to  the  magnitude  of  the  nonlinearity,  significant  error  can  occur 
as  a  result  of  neglecting  the  higher  harmonics  in  the  series  expansion  of  the 
force-displacement  relationship.  The  truncation  of  the  higher  harmonics  is  an 
inherent  drawback  of  the  one-term  describing  function  approach. 

An  experimental  study  of  the  limit  cycle  response  of  aerodynamic  surfaces 
with  structural  nonlinearities  is  discussed  in  Reference  5.  During  this  study 
a  model  was  developed  with  root  structural  behavior  which  gave  a  variety  of 
nonlinear  stiffness  characteristics  to  the  aerosurface.  Wind  tunnel  tests 
were  conducted  and  data  obtained  to  verify  the  model's  response  character¬ 
istics.  It  was  noted  by  the  authors  of  Reference  5  that  there  exists  a  need 
for  improved  analytical  tools  to  accurately  describe  the  nonlinear  behavior  in 
order  to  better  correlate  experimental  and  analytical  results 

The  objective  of  the  present  analytical  study  was  to  develop  a  technique 
to  predict  limit  cycle  response  of  aerosurfaces  with  discrete  structural  non- 
linearities  that  retains  the  flexibility  and  well  defined  procedures  of  the 
describing-function  approach  (References  1  and  2)  yet  provides  greater  ac¬ 
curacy  and  generality  in  modeling  the  nonlinear  system  behavior.  To  meet  this 
objective  an  asymptotic  expansion  method  was  used  to  model  the  nonlinear 
force-displacement  relationship  that  results  when  nonlinearities  of  the  type 
shown  in  Figure  2  are  introduced  at  the  aerosurface  support.  The  primary 
difference  between  the  asymptotic  method  and  the  describing-function  method  is 
the  capability  of  the  asymptotic  method  to  include  higher  harmonics  in  the 
representation  of  the  nonlinearity  and  obtain  successively  higher  order  ap¬ 
proximations  to  the  limit  cycle  response.  Using  the  asymptotic  expansion 
technique  a  broader  category  of  discrete  nonlinearities  than  those  shown  in 
Figure  2  can  be  modeled. 

Specifically,  the  problem  investigated  during  the  present  study  was  the 
limit  cycle  response  of  an  aerodynamic  surface  in  a  subsonic  airstream.  Fig¬ 
ure  3.  The  nonlinearities  shown  in  Figure  2  were  assumed  to  act  at  the  root 
support  springs  and  K  ^  shown  in  Figure  3.  This  problem  is  represenative 
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of  a  missile  control  surface  with  a  loose  hinge  and/or  joint  slippage  in  the 
surface  support  structure  and/or  actuator.  The  aerodynamic  forces  acting  on 
the  surface  were  modeled  using  a  steady  state  aerodynamic  theory.  This  theory 
assumes  the  lifting  force  is  proportional  to  and  in  phase  with  the  torsional 
motion  of  the  surface  which  is  assumed  to  be  sinusoidal.  Simple  aerodynamics 
were  used  so  that  the  influence  of  the  nonlinearities  on  the  surface  response 
could  be  evaluated  by  a  much  more  tractable  computational  effort.  The  effects 
of  using  a  more  sophisticated  aerodynamic  theory  for  the  describing-function 
approach  were  investigated  and  documented  in  Reference  2.  There  it  was  shown 
that  a  more  sophisticated  aerodynamic  theory  can  substantially  change  the 
flutter  results  employed  to  predict  aeroelastic  response,  but  has  no  impact  on 
modeling  the  nonlinear  behavior.  The  authors  believe  that  the  conclusions 
presented  in  Reference  2  are  applicable  to  the  present  analysis  effort  as  the 
linear  aerodynamic  theory  does  not  impact  the  representations  of  the  struc¬ 
tural  nonlinearities  considered  in  this  study. 

This  study  was  organized  into  two  separate  tasks.  The  first  task,  dis¬ 
cussed  in  Section  2.0,  involved  development  of  the  asymptotic  approximations 
for  the  limit  cycle  response  for  each  of  the  two  nonlinear ities  shown  in  Fig¬ 
ure  2.  First  and  second  order  asymptotic  solutions  were  formulated  and  com¬ 
pared  with  the  results  obtained  using  the  describing-function  method  of  Re¬ 
ferences  1  and  2.  The  equations  for  including  the  higher  order  harmonics  in 
the  time  history  response  calculations  were  also  derived  under  this  task.  The 
second  task,  discussed  in  Section  3,  involved  application  of  the  developed 
asymptotic  solutions  to  predict  the  limit  cycle  response  of  the  baseline  aero¬ 
dynamic  surface  chosen  for  this  study.  This  baseline  surface  design  is  based 
on  the  Harpoon  missile  control  surface  and  is  shown  in  Appendix  A.  In  Appen¬ 
dix  A,  representative  flutter  results  for  selected  root  spring  configurations 
are  presented.  These  results  are  used  in  the  applications  discussed  in  Sec¬ 
tion  3.  Appendix  8  contains  the  detailed  derivation  of  the  asymptotic  solu¬ 
tions.  While  key  equations  are  presented  in  the  text  of  Section  2,  inter¬ 
mediate  steps  in  the  derivations  are  detailed  in  Appendix  B. 
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2.0  DEVELOPMENT  OF  ASYMPTOTIC  SOLUTION 


Application  of  the  method  of  asymptotic  expansions  to  model  nonlinear 
behavior  of  an  aerodynamic  surface  response  forms  the  basis  of  the  procedures 
developed  in  this  investigation.  The  asymptotic  expansion  method  used  is  the 
Krylov-Bogoliubov-Mitropolski  (KBM)  technique.  This  technique  is  a  perturbation 
method  based  on  a  more  general  approach  known  as  the  method  of  averaging,  and 
is  discussed  in  detail  in  References  6  and  7.  In  the  method  of  averaging,  the 
motion  is  assumed  to  vary  slowly  with  time  and  the  amplitude  and  phasing  are 
determined  as  time  dependent  functions  of  the  system  nonlinearities.  The 
advantage  of  the  method  of  averaging,  and  in  particular  the  KBM  technique,  is 
the  stability  of  the  solutions  over  sufficiently  long  intervals  of  time.  This 
type  of  solution  is  necessary  for  quantitative  as  well  as  qualitative  analysis 
of  the  limit  cycle  behavior  of  a  nonlinear  system.  Another  advantage  of  the 
asymptotic  expansion  technique  is  that  it  provides  a  more  accurate  determination 
of  the  nonlinear  load-displacement  relationship  than  the  one-term  describing- 
function  approach  of  References  1  and  2,  yet  can  be  used  in  much  the  same 
manner.  In  both  the  asymptotic  and  describing-function  techniques,  the  derived 
expressions  for  the  load-displacement  relationship  are  used  to  define  an 
"effective"  stiffness  for  the  nonlinear  element.  This  effective  stiffness  is 
used  in  subsequent  aeroelastic  analyses. 

The  procedures  for  investigating  the  limit  cycle  response  developed  in  this 
study  require  defining  the  flutter  characteristics  of  the  aerodynamic  surface. 

In  this  study,  flutter  analysis  of  both  rigid  and  flexible  configurations  were 
analyzed.  The  flutter  analysis  is  performed  using  the  linearized  form  of  the 
equations  of  motion,  as  was  done  in  the  describing-function  and  other  methods. 
References  1  through  4.  This  approach  assumes  that  an  effective  stiffness 
adequately  describes  the  elastic  behavior  for  the  infinitesimal  oscillations  at 
flutter  onset.  The  details  of  the  flutter  analysis  for  the  baseline  surface 
analyzed  in  this  study  are  presented  In  Appendix  A. 

During  this  study,  the  nonlinear  aeroelastic  system  was  represented  by  the 
second  order  differential  equation 

2 

^-4  +  w  2x  -  c  f  (x)  +  Q  (x)  (1) 

dt2  0 


7 


In  Equation  (l),e  f(X)  represents  the  nonlinear  force  acting  on  the  elastic 
system,  X  is  the  surface  root  displacement,  <*»0  is  the  natural  frequency  of  the 
linear  system  («a0)  and  Q(x)  is  the  aerodynamic  forcing  function.  The 
f (X)  term  appears  in  Equation  (1)  as  a  result  of  the  nonlinear  load-displace¬ 
ment  relationship  for  the  root  support  springs  of  the  surface.  It  is  the 
nonlinear  force  term  that  is  used  to  determine  the  form  of  the  asymptotic 
expansions  derived  to  approximate  the  limit  cycle  response  of  the  nonlinear 
system. 


In  the  asymptotic  method  a  perturbation  technique  is  used  to  expand  the 
contribution  of  the  nonlinearity  in  terms  of  a  small  parameter  or  gage  func¬ 
tion,*  .  This  expansion  takes  the  form  of  an  asymptotic  series  comprised  of 
integer  powers  of  «  .  For  the  system  respresented  by  Equation  (1)  the  form  of 
the  asymptotic  solution  defined  by  the  KBM  technique  is  given  by 

x  *  A  cos  *  +  £  en  U  (A,i|>)  +  0  (eN+1)  (2) 

n*l  n 


For  the  freeplay  nonlinearity  this  equation  was  used  directly.  For  the  preload 
nonlinearity  it  was  modified  as  will  be  discussed  in  Section  2.2.  The 
displacement,  X,  as  defined  in  Equation  (2)  consists  of  a  linearly  independent 
combination  of  periodic  functions,  the  first  term  being  the  fundamental  har¬ 
monic.  The  remaining  terms  expanded  in  powers  cf  e  represent  the  asymptotic 
approximation  of  the  nonlinear  contribution  to  the  response.  The  functions  Un 
are  periodic  functions  comprised  of  higher  harmonics  of  the  phasing  parameter  ^ 
and  the  surface  root  amplitude  A.  N  indicates  the  order  of  the  asymptotic 
approximation  and  the  symbol  0(*N+1)  represents  term  of  order  greater  than  N. 

A  and  are,  in  general,  functions  of  time  defined  by  the  ordinary  differen¬ 


tial  equations 


4.  t 


*n  +  °  <'N+1) 


d<ii 

7F 


.  v'  n  „  .  n  [  N+l * 
“o  +  L  e  6  +  0  (e  ) 
0  n«l  n 


(4) 


The  right  hand  side  of  Equations  (3)  and  (4)  are  series  expansions  of  €,  the 
gage  function,  and  parameters  an  and/3n.  These  an  and /3n  terms  are  determined 
by  integrating  the  nonlinear  force  term  f ( X )  that  appears  in  Equation  (1)  over 
the  total  period  of  the  oscillation  which  is  defined  to  be  27r  .  The  amplitude 
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and  phasing  are  thereby  time  dependent  functions  of  the  nonlinear  load- 
displacement  relationship  defined  by  Equations  (3)  and  (4). 

In  the  asymptotic  method,  as  in  most  perturbation  methods,  the  linear 
system  results  are  approached  as  «  approaches  zero.  In  this  case,  the  ampli¬ 
tude,  A,  becomes  constant  and  the  phasing,  <|» ,  becomes  a  linear  function  of 
time  with  a  constant  frequency,  Equations  (3)  and  (4).  The  displacement  in 
Equation  (2)  then  takes  the  form  of  simple  harmonic  motion  of  linear  elastic 
systems.  The  basis  of  the  asymptotic  solution  approach  is  to  determine  the 
solutions  of  equations  (2)  through  (4)  for  successively  higher  orders  of  <• 

The  procedure  for  determining  the  asymptotic  solutions  depend  on  deriving 
appropriate  functional  forms  of  Un,  an,  and  /3n.  This  is  accomplished  by  ex¬ 
pressing  the  load  developed  in  the  nonlinear  root  springs  by 

F  (x)  *  K(x)X  (5) 

K( X)  is  the  nonlinear  stiffness  associated  with  the  root  support  spring  and  X 
is  the  root  displacement  defined  by  Equation  (2).  To  determine  the  Un,  an  and 
coefficients  only  the  first  term  of  Equation  (5)  is  retained  and  this  term 
is  expanded  in  the  Fourier  Series.  The  functions  «n,  /3n  and  Un  are  expressed 
in  terms  of  the  coefficients  of  this  Fourier  Series  and  therefore  are  inte¬ 
grals  of  the  nonlinear  force  term  evaluated  over  a  period  of  2 ir.  In  this 
manner,  the  equations  defining  the  asymptotic  solutions  of  Equations  (2) 
through  (4)  are  completely  defined.  During  this  study  the  asymptotic  solu¬ 
tions  were  obtained  for  first  and  second  order  approximations,  i.e.,  terms  in 
Equations  (3)  and  (4)  up  to  «2.  The  details  of  the  Fourier  Series  expansions 
and  the  determination  of  the  asymptotic  solutions  are  given  in  Appendix  B. 

The  results  of  the  asymptotic  expansion  for  the  two  types  of  nonlinearities 
considered  in  this  study  are  summarized  below. 

2.1  Freeplay  Nonlinearity 

The  freeplay  nonlinearity  is  illustrated  in  Figure  2(a).  The  asymptotic 
solution  for  the  freeplay  nonlinearity  is  derived  directly  from  the  equations 
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developed  in  the  preceding  section.  The  waveform  of  the  developed  load  will 
take  one  of  the  two  shapes  shown  in  Figure  4,  depending  on  the  relationship 
between  the  magnitudes  of  the  freeplay,  S,  and  the  amplitude.  A,  of  displace¬ 
ment.  For  A  less  than  S  no  load  is  developed.  Figure  4(a).  When  A  is  greater 
than  S  the  load  is  as  shown  in  Figure  4(b).  The  nonlinear  load  is  then  de¬ 
termined  by  specifying  the  freeplay  S  and  the  amplitude  A.  The  approach  used 
throughout  this  study  is  based  on  derivation  of  the  "effective"  stiffness  K  of 
the  nonlinear  spring.  This  effective  stiffness  is  used  in  Equation  (5)  to 
determine  the  corresponding  load.  With  the  effective  load  substituted  in 
Equation  (1),  the  equivalent  linearized  system  is  analyzed  by  conventional 
methods.  It  is  the  determination  of  this  effective  stiffness  term  that  uti¬ 
lizes  the  asymptotic  expansion  techniques. 

The  "effective"  stiffness  K  of  the  nonlinear  spring  including  the  in¬ 
fluence  of  freeplay  is  defined  in  the  asymptotic  method  as: 

IT  *  K  ^  1  +  e  B2  +  •  • .  B^  (6) 

In  this  expression,  the  Bn  terms  are  directly  related  to  the  Pn  coefficients 
of  the  expansion  of  Equation  (4).  The  relationship  between  the/3n  and  Bn 
terms  is  derived  in  Appendix  B.  Therefore,  rather  than  actually  solving  Equa¬ 
tion  (4),  once  the  coefficients  are  determined  and  therefore  the  Bn's,  the 
"effective"  stiffness  is  computed  from  Equation  (6)  directly.  For  the  free¬ 
play  nonlinearity  all  on  coefficients  are  zero  due  to  the  symmetry  of  the 
freeplay  nonlinearity  and  the  single  valued  odd  form  of  the  nonlinear  load. 

This  implies,  from  Equation  (3),  that  the  amplitude  of  the  limit  cycle  oscillation 
is  constant. 

The  expression  for  the  first  order  approximation  of  the  effective  stiff¬ 
ness,  hereafter  referred  to  as  the  first  order  solution,  contains  only  the  B] 
term  in  Equation  (6).  The  form  of  this  term  (from  Appendix  B)  is  given  as 

e  B1  ■  0  (A<S) 

6  Bi  ■  ir  0  -  2  *1  +  cos  tl  sin  2  (7) 

t.j  »  cos"1  (S/A)  (8) 
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Using  Equations  (6),  (7)  and  (8)  the  relationship  between  the  effective  stiff¬ 
ness  and  the  linear  spring  rate  for  a  freeplay  nonlinearity  may  be  obtained  to 
the  first  order. 

The  form  of  Bi  developed  from  the  asymptotic  expansion  technique  is  iden¬ 
tical  to  that  developed  using  the  describing-function.  References  1  and  2. 

Since  the  describing  function  results,  and  hence  the  first  order  solutions  are 
discussed  in  detail  in  References  1  through  4,  they  will  only  be  summarized 
here  for  comparison  and  the  remaining  discussion  focused  on  the  higher  order 
solutions. 

The  second  order  solution  is  obtained  by  adding  the  effect  of  second  order 
terms  to  the  first  order  solution.  In  order  to  determine  the  second  order 
approximation  of  the  effective  stiffness,  it  is  necessary  to  derive  the  ex¬ 
pression  for  the  82  coefficient  that  appears  in  Equation  (6).  Using  the  pro¬ 
cedures  for  the  asymptotic  method,  the  second  order  correction  term  for  a 
single  degree  of  freedom  with  a  freeplay  nonlinearity  is  of  the  form  (from 
Appendix  B) 


e2  B2  -  0  (A<S) 

-2  O  3.1  V  -Is-  l"r  -  -  cos  t.  sin  nt^ 
'  “2  *  n-3,5  77  L  "  "  1 


(9) 


where 


sin  (n-l)t^  sin  (n+l)t.j' 
'  2  (n-l  F  +  2  (n+l )  . 


(10) 


1? 


and  t)  is  definded  In  Equation  (8).  The  second  order  correction  terms  of 
Equation  (9)  contains  the  coefficients  of  the  fundamental  harmonic  in  the 
Fourier  Series  expansion  of  the  load  and  coefficients  of  all  higher  order 
harmonics.  Therefore  the  second  order  solutions  are  not  restricted,  as  is  the 
describing-function  approach  of  References  1  through  4,  to  a  one-term,  first 
harmonic  approximation.  Substituting  Equations  (7)  through  (10)  into  Equa¬ 
tion  (6),  the  relationship  between  the  effective  stiffness  and  the  linear 
spring  rate  for  a  freeplay  nonlinearity  may  be  obtained  for  the  second  order 
approximation.  This  relationship  is  shown  in  Figure  5  as  a  function  of  the 
amplitude  of  motion  to  freeplay  ratio,  (A/S).  For  amplitude  ratios  (A/S)  less 
than  1,  the  effective  stiffness  is  zero.  As  the  amplitude  increases,  the 
magnitude  of  K  approaches  that  of  linear  stiffness  K.  This  corresponds  physi¬ 
cally  to  the  nonlinearity  becoming  less  and  less  significant  as  the  (A/S) 
ratio  increases  and  «  approaching  zero  in  Equations  (3)  and  (4). 

As  shown  in  Figure  5,  there  is  little  difference  in  the  results  of  the 
first  and  second  order  solutions  where  the  nonlinearity  is  most  significant, 
low  (A/S)  values,  and  the  solutions  converge  very  rapidly  as  (A/S)  increases. 
The  implication  here  is  that  the  stiffness  behavior  of  system  with  a  freeplay 
nonlinearity  is  dominated  almost  completely  by  the  fundamental  harmonic  and 
the  second  order  corrections  have  little  impact  on  the  response.  For  a  single 
degree  of  freedom  system,  the  effective  stiffness  value  can  be  expressed  in 
terms  of  an  effective  frequency  of  the  nonlinear  system.  For  the  stiffness 
ratios  shown  in  Figure  5,  the  corresponding  first  and  second  order  solutions 
expressed  in  terms  of  frequency  are  shown  in  Figure  6.  In  this  figure  the 
asymptotic  solutions  are  compared  with  numerical  simulation  results  obtained 
for  a  single  degree  of  freedom  system  in  Reference  1.  As  the  figure  indicates, 
the  asymptotic  solutions  predict  frequency  ratios  very  close  to  those  given  by 
the  numerical  simulation. 

2.2  Preload  Nonlinearity 

For  the  preload  nonlinearity  shown  in  Figure  2(b),  the  asymptotic  func¬ 
tions  or  displacements,  defined  by  Equation  2,  were  modified  to  account  for  the 
non-symmetry  of  the  load  displacement  relationship.  This  displacement  func- 
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FIGURE  5  EFFECTIVE  STIFFNESS  RATIO  FOR  A  FREEPLAY  NONLINEARITY 


FIGURE  6  COMPARISON  OF  ASYMPTOTIC  AND  NUMERICAL  SIMULATION  RESULTS 
FOR  A  FREEPLAY  NONLINEARITY 


tion  was  assumed  to  be  of  a  similar  form  to  that  used  in  Reference  1  and  is 


A  +  A,  cos  ifr  +  E  un  +  0  (e  ^ 

o'  n-1 


OD 


The  coefficients  A0  and  Ai  were  defined  such  that  the  energy  stored  in  the 
nonlinear  spring  is  the  same  for  both  positive  and  negative  displacements, 
addition,  it  was  required  that  Equation  (11)  result  in  a  positive  amplitude 
equal  to  the  initial  displacement.  Thus,  the  amplitude  coefficients  that 
appear  in  Equation  (11)  are  obtained  from 


A,  «  |  +  \  ^2  PA 


(P<A<P+2S) 


A1  *7 


(A  -  2S)Z  +  4PS  (A>P+2S) 

In  both  cases,  the  coefficient  A0  is  obtained  from  the  relationship 


(12) 

(13) 


In 


A.  +  A 
1  o 


(14) 


It  has  been  assumed  that  the  influence  of  a  preload  nonlinearity  is  related  to 
positive  displacements  of  the  system.  When  amplitude  of  motion.  A,  is  la;s 
than  the  preload  P,  A]  equals  A  and  A0  is  zero.  For  this  situation,  the  non¬ 
linear  problem  is  reduced  to  a  linear  problem. 


The  waveform  of  the  developed  load  in  the  nonlinear  spring  will  take  the 
shapes  shown  in  Figure  7.  As  before,  these  waveforms  are  dependent  on  the 
relationship  between  the  magnitudes  of  the  freeplay,  preload  and  amplitude  of 
motion.  Proceeding  as  was  done  for  the  freeplay  case,  the  coefficients  of 
Equation  (6)  were  defined  using  the  load  displacement  relationship  illustrated 
in  Figure  7.  The  first  order  approximation  to  the  "effective"  stiffness  is  of 
the  form 

c  B1  *  1.0  (A<P)  (15) 

*  8i  ■  t  [’  -  +  h;  <p-V  s1n  *i  - 1 s1" 2  4i]  (,6) 

(P<A<P+2S) 

For  amplitude  of  motion  A>P+2S  we  have, 

c  8,  .i[.  +  t,  -  t2  -  ^  (P*2S-V  sin  t,  *  (17) 

(P-A,)  sin  t2  *  {  (sin  2  t,  -  sin  2  t2)] 
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where 


(18) 

09) 


As  was  determined  for  the  freeplay  case  the  first  order  asymptotic  solution 
and  the  describing-function  solution  of  References  1  and  2  are  equivalent. 
Proceeding  to  the  second  order  approximation  the  B2  coefficient  of  Equation  (6) 
is  of  the  form  (from  Appendix  B) 


e2  B, 


,  «  1 

_  1  r  n _ _ 

n  n=2  n^T 


A  ) 
o' 


sin  nt 


1 


r-J! 


(20) 


rn  is  defined  by  Equation  (10),  t]  by  Equation  (16)  and  B]  is 
tion  (15).  For  the  amplitudes  A>P+2S 


2 

c 


,  ** 


n-2  n  -1 


given  by  Equa 


(21) 


-^2  cos  t-j  sin  nt^  +  on 
B]  and  t2  are  defined  by  Equations  (17)  and  (19)  and 


T  sin  (n-l)t£  -  sin  (n-l)t-|  sin  (n+l)t^  -  sin  (n+ljt^ 
an  “  L  2"  (n-1 )  +  2  (n+1 )  J  ^ 

The  first  and  second  order  solutions  for  the  stiffness  ratio  is  plotted 
in  Figure  8  as  a  function  of  the  amplitude  of  motion  to  freeplay  ratio  for  a 
freeplay  to  preload  ratio  (S/P)  of  one.  For  amplitudes  of  motion  less  than 
the  preload  P  the  frequency  ratio  is  one  and  the  response  is  linear.  As  the 
amplitudes  of  motion  increases,  the  stiffness,  and  in  turn,  the  frequency 
decreases.  This  softening  response  is  due  to  the  deadspace  in  the  spring 
causing  the  effective  stiffness  to  be  less  than  the  linear  value.  As  the 
amplitude  increases  well  beyond  the  nonlinear  region,  the  influence  of  the 
nonlinearity  becomes  small  and  the  magnitude  of  frequency  coefficient  ap¬ 
proaches  one. 
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AMPLITUDE 

FREEPLAY 


FIGURE  8  EFFECTIVE  STIFFNESS  RATIO  FOR  A  PRELOAD  NONLINEARITY 


While  both  first  and  second  order  solutions  shown  in  Figure  8  exhibit  the 
same  trends,  the  second  order  solution  predicts  a  significantly  greater  re¬ 
duction  in  stiffness  for  amplitudes  of  motion  near  P+2S  where  the  nonlinearity 
is  most  significant.  At  higher  values  of  the  amplitude,  the  first  and  second 
order  solutions  converge  and  asymptotically  approach  the  linear  solution.  The 
difference  in  the  two  solutions  near  the  nonlinear  region  is  directly  a  result 
of  including  higher  harmonics  in  the  second  order  approximation.  These  har¬ 
monics  are  seen  to  influence  the  response  most  when  the  amplitudes  of  motion 
fall  within  the  deadband  region  of  the  load  displacement  curve.  Figure  7. 

A  comparison  between  the  asymptotic  solutions  and  numerical  solutions  are 
shown  in  Figure  9.  The  numerical  solutions.  Reference  1,  were  obtained  by 
directly  integrating  the  single  degree  of  freedom  equations  of  motion.  Fi¬ 
gure  9  shows  that  the  first  order  approximation  (describing-function)  is  un¬ 
conservative  in  predicting  the  effective  frequency  behavior  when  compared  to 
the  numerical  results  for  amplitudes  near  the  preload-plus-freeplay  region, 
(A/S=3).  On  the  other  hand,  the  second  order  solution  tends  to  be  conserva¬ 
tive  in  predicting  the  effective  stiffness.  At  amplitudes  of  motion  near  or 
greater  than  twice  the  preload-plus-freplay  values,  there  is  little  difference 
between  the  numeri:al  and  first  or  second  order  solutions.  Based  on  these 
comparisons,  it  is  apparent  that  the  higher  harmonics  contribute  significantly 
to  the  response  at  amplitudes  of  motion  near  the  nonlinear  region  of  the 
load-displacement  relationship  for  the  preload  nonlinearity. 

2.3  Time  History  Results 

In  addition  to  the  computation  of  the  effective  stiffness  of  the  non¬ 
linear  system,  the  asymptotic  method  provides  a  means  of  obtaining  higher 
order  approximations  to  the  response  time  history.  The  time  dependent  motion 
of  the  nonlinear  system  as  defined  in  the  asymptotic  method  is  given  by  Equa¬ 
tion  (2).  In  this  expression,  the  functions  Un  contain  the  contributions  of 
higher  harmonics  in  the  response  and  exclude  any  contribution  of  the  funda¬ 
mental  harmonic.  Application  of  the  asymptotic  method  consists  of  determining 
the  appropriate  expression  for  Un  based  on  the  order  of  the  asymptotic  ap¬ 
proximation  and  the  number  of  harmonics  desired.  The  general  form  for  the  ex¬ 
pression  for  each  type  of  nonlinearity  considered  in  this  study  is  given  be¬ 
low. 
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FIGURE  9  COMPARISON  OF  ASYMPTOTIC  AND  NUMERICAL  SIMULATION  RESULTS 
FOR  A  PRELOAD  NONLINEARITY 
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2.3.1  Freeplay  Nonlinearity 

For  the  freeplay  nonlinearity  the  form  of  the  U  function  is  determined 
using  the  coefficients  of  the  Fourier  Series  expansion  of  the  load  displace¬ 
ment  relationship.  From  the  results  of  Appendix  B  to  the  second  order 

U,  »  -  —  £  r„  -  \  cos  t,  sin  ntJ  cos  nut  (23) 

1  *  n=3,5  L  n  n  1  II 

N  is  odd,  Tn  is  defined  in  Equation  (10)  and  w  is  the  effective  frequency. 

The  time  history  results  for  the  freeplay  nonlinearity  are  shown  in  Fig¬ 
ures  10  and  11.  In  Figure  10,  the  time  history  response,  defined  by  Equa¬ 
tion  (2),  is  plotted  for  an  A/S  ratio  of  3  and  an  initial  amplitude  of  0.2°. 
The  limit  cycle  motion  shown  in  Figure  10  exhibits  a  very  regular  periodic 
motion  with  an  almost  constant  amplitude  and  frequency  corresponding  to  effec¬ 
tive  stiffness  value  determined  by  the  relationship  defined  in  Paragraph  2.1. 
The  period  of  the  time  history  indicates  that  the  response  is  controlled  al¬ 
most  entirely  by  the  fundamental  harmonic,  the  first  term  of  Equation  (2). 

The  time  history  of  the  second  term  of  Equation  (2),  corresponding  to  the 
contribution  of  higher  harmonics,  is  shown  in  Figure  11.  Terms  up  to  the 
eleventh  harmonic  were  included  in  Equation  (23)  in  the  determination  of  U. 
Note  that  although  the  frequency  of  U  is  higher,  the  amplitude  is  signifi¬ 
cantly  lower  than  that  of  fundamental  harmonic.  These  results  are  consistent 
with  the  comparisons  between  the  first  and  second  order  solutions  discussed  in 
Paragraph  2.1.  The  results  show  the  first  order  solution  utilizing  a  one  term 
Fourier  Series  expansion  of  the  load  is  sufficient  for  the  freeplay  nonlin¬ 
earity  and  the  higher  harmonics  of  the  second  order  solutions  have  little 
impact  on  the  system  response. 


2.3.2  Preload  Nonlinearity 


For  the  preload  nonlinearity  the  form  of  the  U  function  is  determined  in 
the  same  manner  as  was  done  for  the  freeplay  case.  From  the  results  of  Ap¬ 
pendix  B  the  second  order  solutions  for  amplitudes  of  motion  in  the  range  P<A< 


P+2S  is  defined  as 
ui  ■ 


2fcTi  *  (2"t') 


-  2  sin 


■  jz  -y 

’  n-2  n  -1 


cos  nwt 


in  t,l 

h  [*’’■*«> s,n  "*1 


‘"I! 


2  2 


(24) 
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FIGURE  10  LIMIT  CYCLE  TIME  HISTORY:  FUNDAMENTAL  PLUS  HIGHER  HARMONICS 
FREEPLAY  NONLINEARITY 
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FIGURE  11  TIME  HISTORY:  HIGHER  HARMONIC  -  FREEPIAY  NONLINEARITY 
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rn  and  t]  are  given  in  Equations  (10)  and  (16).  For  amplitudes  of  motion  A> 
P+2S 


7  [£?  (‘i +  ’  -  '*) +  (s1"  ' s,n  **) +  (V‘i) 

-  sin  nt^j- 


(25) 


2  cos  t^  sin  nt^  +  on| 


1 1 ,  t2  and  crn  are  defined  by  Equations  (18),  (19)  and  (22). 


The  time  history  results  for  the  preload  cases  are  shown  in  Figures  12 
and  13.  In  Figure  12  the  displacement  time  history  for  a  preload  nonlinearity 
with  an  initial  amplitude  of  0.2°  is  shown.  In  this  figure  the  two  distin¬ 
guishing  features  of  the  preload  time  history  is  the  offset  of  the  response 
from  zero  and  the  presence  of  higher  harmonics  in  the  response.  The  offset  is 
due  to  the  presence  of  the  preload  and  the  unsymmetric  nature  of  the  nonlinear 
load-displacement  relationship.  The  higher  harmonics  in  the  amplitude  are  due 
to  the  contribution  of  terms  appearing  in  1).  The  time  history  plot  of  the 
higher  harmonics  in  the  U  coefficient  is  shown  in  Figure  13.  In  Figure  13  the 
offset  in  the  response  is  also  present.  Also  note  that  comparing  the  magni¬ 
tude  of  the  time  histories  in  Figures  2  and  13  shows  that  the  U  terms  are 
greater  than  forty  percent  of  the  total  amplitude.  This  indicates  that  a 
significant  amount  of  the  strain  energy  developed  in  the  root  support  springs 
cycles  with  a  frequency  greater  than  the  frequency  of  the  fundamental  harmon¬ 
ic.  Based  on  these  observations,  the  influence  of  higher  harmonics  is  seen  to 
be  significant  in  the  prediction  of  both  the  effective  frequency  and  limit 
cycle  response  of  the  system  response  for  a  preload  nonlinearity. 
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FIGURE  12  LIMIT  CYCLE  TIME  HISTORY:  FUNDAMENTAL  PLUS  HIGHER  HARMONICS 
PRELOAD  NONLINEARITY 
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3.0  DETERMINATION  OF  AERODYNAMIC  SURFACE 
LIMIT  CYCLE  RESPONSE 


In  this  section  the  equations  and  techniques  developed  in  Section  2.0  for 
the  freeplay  and  preload  nonlinearities  are  applied  to  predict  the  limit  cycle 
response- of  the  baseline  aerodynamic  surface.  The  frequency  and  dynamic  pres¬ 
sure  at  which  the  limit  cycle  oscillations  are  sustained  were  determined  as  a 
function  of  the  amplitude  of  motion  and  magnitude  of  the  nonlinearity.  The 
aeroelastic  data  used  was  based  on  linear  flutter  results  obtained  for  the 
effective  system  where  the  nonlinear  terms  in  the  equations  of  motion  were 
replaced  by  the  corresponding  effective  stiffness  calculated  by  the  asymptotic 
methods.  The  flutter  results  for  both  a  rigid  and  flexible  representation  of 
the  baseline  aerodynamic  surface  were  determined  using  standard  eigenanalysis 
procedures  and  are  presented  in  Appendix  A. 

The  geometry  and  physical  characteristics  of  the  baseline  control  surface 
are  also  presented  in  Appendix  A.  The  results  of  this  section  and  Appendix  A 
are  for  a  specific  aerodynamic  surface  using  a  simplified  aerodynamic  theory. 
However,  the  application  of  the  developed  procedure  is  not  restricted  to  these 
conditions.  The  procedures  are  applicable  to  a  variety  of  surface  geometries, 
aerodynamic  theories  or  flutter  analysis  techniques. 

Application  of  the  asymptotic  expansion  results  to  predict  the  limit  cycle 
response  of  the  baseline  aerodynamic  surface  follows  the  same  procedure  as 
developed  in  References  1  and  2.  The  effective  uncoupled  frequency  was 
calculated  from  the  asymptotic  solutions  for  each  nonlinearity  and  corresponding 
degree  of  freedom.  The  dynamic  pressure  defined  by  Equation  26,  at  which  the 
limit  cycle  oscillation  is  sustained,  was  then  determined  from  the  data  in 
Appendix  A. 


0  *  \  p  U2  (26) 

The  results  for  the  freeplay  and  preload  nonlinearities  considered  in  this  study 
follow. 
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3.1  FREEPLAY  NONLINEARITY 


The  procedure  followed  to  determine  the  limit  cycle  response  of  the  aero¬ 
dynamic  surface  with  root  freeplay  structural  nonlinearities  is  shown  in  Fig¬ 
ure  14.  This  figure  shows  the  step-by-step  procedure  for  introducing  the 
nonlinear  effects  in  either  or  both  root  pitch  and  root  roll  degrees  of  free¬ 
dom. 


Results  for  the  rigid  baseline  control  surface  with  a  freeplay  nonlin¬ 
earity  in  the  root  pitch  degree  of  freedom  are  shown  in  Figure  15.  These 
results  are  for  a  freeplay  nonlinearity  in  the  root  pitch  degree  of  freedom, 

S.,  of  0.2  degrees  and  an  uncoupled  root  pitch  frequency  of  21 5H2 .  The  data 

<D 

indicates  the  variation  in  the  effective  dynamic  pressure  at  which  limit  cycle 
oscillations  will  be  sustained  as  a  function  of  the  ampl itude-to -freeplay 
ratios  in  the  root  pitch  degree  of  freedom.  This  figure  shows  the  influence 
of  the  freeplay  nonlinearity  is  most  pronounced  for  amplitudes  of  motion  near 
the  freeplay  value.  For  successively  larger  amplitudes  of  motion,  the  effec¬ 
tive  stiffness  approaches  the  uncoupled  stiffness  value  and  the  critical  dy¬ 
namic  pressure  approaches  that  of  the  linear  system. 

The  results  obtained  using  the  first  and  second  order  asymptotic  solu¬ 
tions  to  predict  the  limit  cycle  behavior  are  compared  in  Figure  15.  They 
indicate,  as  was  found  in  Section  2.0,  that  the  first  and  second  order  solu¬ 
tions  differ  only  slightly  for  the  freeplay  case.  Again,  this  behavior  is 
attributed  to  the  dominance  of  the  fundamental  harmonic  in  the  response. 

Results  for  freeplay  nonlinearities  in  root  roll  and  for  two  freeplay 
nonlinearities  were  found  to  differ  only  slightly  from  those  of  References  1  and  2. 
The  results  for  these  cases  are  not  repeated  here  and  the  reader  is  referred 
to  References  1  and  2  for  more  details  on  the  influence  of  freeplay 
nonlfnearities  on  the  limit  cycle  behavior. 
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FIGURE  15  DYNAMIC  PRESSURE  TO  SUSTAIN  LIMIT  CYCLE  OSCILLATION;  RIGID 
SURFACE  WITH  ROOT  PITCH  FREEPLAY  NONLINEARITY 
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3.2  PRELOAD  NONLINEARITY 


The  procedure  used  for  predicting  the  limit  cycle  response  of  the  rigid 
baseline  control  surface  with  a  preload  nonlinearity  is  very  similar  to  that 
presented  in  the  preceding  section  for  the  freeplay  nonlinearity.  The  compu¬ 
tational  steps  to  be  followed  to  predict  the  limit  cycle  response  are  shown  in 
Figure  16.  Linear  system  flutter  analyses  are  conducted  for  variations  in  the 
effective  stiffness  parameter.  These  flutter  results  are  then  modified  to 
account  for  the  presence  of  the  structural  nonlinearities. 

Using  the  second  order  asymptotic  solutions,  the  results  shown  in  Fig¬ 
ure  17  were  obtained  for  a  rigid  control  surface  having  a  single  root  roll 
preload  nonlinearity.  This  figure  shows  the  change  in  the  dynamic  pressure  as 
a  function  of  root  roll  amplitude  of  motion  to  freeplay  ratio,  A/S,  for  vary¬ 
ing  freeplay-to-preload,  S/P,  ratios.  For  amplitudes  of  motion  less  than  the 
preload,  the  critical  dynamic  pressure  equals  the  flutter  dynamic  pressure  of 
the  linear  system.  As  the  amplitude  of  motion  increases,  the  influence  of  the 
freeplay  is  reflected  in  the  rise  of  the  dynamic  pressure.  This  occurs  as  a 
result  of  the  softening  effect  on  the  effective  root  roll  stiffness  which 
results  in  a  higher  dynamic  pressure  for  this  particular  aerodynamic  surface. 
As  the  amplitude  of  motion  continues  to  increase,  the  influence  of  the  non¬ 
linearity  decreases  and  the  results  again  approach  those  of  the  linear  system. 

In  Figure  18  the  results  obtained  for  both  first  and  second  order  asymp¬ 
totic  solutions  are  shown  for  a  preload  nonlinearity  in  the  root  pitch  degree 
of  freedom  for  a  rigid  control  surface.  From  this  figure  the  influence  of  the 
higher  harmonics  is  apparent.  While  the  trends  are  basically  the  same,  the 
second  order  solution  predicts  a  lower  value  of  critical  dynamic  pressure  to 
sustain  the  limit  cycle  oscillation  than  does  the  first  order  solution  for 
amplitudes  of  motion  near  the  P+2S  value.  This  is  directly  related  to  the 
effective  stiffness  behavior  predicted  for  this  amplitude  range  as  was  dis¬ 
cussed  in  Section  2.0  and  illustrated  in  Figure  8. 

First  and  second  order  solutions  for  the  rigid  baseline  control  surface 
having  preload  nonlinearities  in  both  root  degrees  of  freedom  are  shown  in 
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FIGURE  16  COMPUTATIONAL  PROCEDURE  FOR  A  PRELOAD  NONLINEARITY 


Figure  19.  The  results  presented  are  for  a  0.2  degree  freeplay  region  in  both 
degrees  of  freedom  and  a  freeplay-to-preload, S/P,  ratio  of  two.  The  two  pair 
of  curves  shown  in  this  figure  are  for  two  different  values  of  the  amplitude 
to-freeplay  ratio,  S/A.  The  effective  stiffness  of  the  preload  nonlinearity 
is  a  double-valued  function  as  illustrated  in  Figure  8.  Therefore,  the  re¬ 
sults  shown  in  Figure  19  are  for  double-valued  amplitude  of  motion  ratios. 

The  larger  S/A  values  correspond  to  amplitudes  less  than  the  quantity  P+2S, 
whereas  the  lower  ratios  correspond  to  amplitudes  in  excess  of  this  value. 

As  for  the  case  of  a  single  nonlinearity,  the  presence  of  the  higher  har¬ 
monics  can  substantially  change  the  prediction  of  the  limit  cycle  response. 

In  both  Figures  18  and  19,  the  nonlinearity  in  root  pitch  is  more  critical  in 
terms  of  initiating  a  limit  cycle  response.  For  root  pitch  nonl inearities  the 
results  indicate  that  the  limit  cycle  response  amplitude  can  be  several  times 
greater  than  the  magnitude  of  the  freeplay  and  can  be  sustained  at  dynamic 
pressures  well  below  the  flutter  critical  value. 

A  flexible  control  surface  having  a  preload  nonlinearities  in  both  root 
degrees  of  freedom  was  also  studied.  The  results  for  a  0.2  degree  freeplay  in 
both  root  degrees  of  freedom  and  a  freeplay-to-preload  ratio, S/P,  of  two  are 
shown  in  Figure  20.  As  with  the  rigid  control  surface  studies  of  Figures  18 
and  19,  the  amplitude  ratios  are  double  valued. 

Figure  20  shows  that  at  higher  uncoupled  roll  frequencies,  the  second 
order  solutions  predict  a  somewhat  higher  critical  dynamic  pressure  than  does 
the  first  order  solution  at  the  same  ampl itude -to -freeplay  ratio.  This  be¬ 
havior  is  a  direct  result  of  the  linear  flutter  data,  as  the  trends  in  the 
effective  stiffness  are  the  same  for  both  the  rigid  and  flexible  cases. 


(ISd) 


FIGURE  19  DYNAMIC  PRESSURE  TO  SUSTAIN  LIMIT  CYCLE  OSCILLATION;  RIGID 
SURFACE  WITH  PRELOAD  NONLINEARITIES  IN  ROOT  PITCH  AND  ROLL 
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FIGURE  20  DYNAMIC  PRESSURE  TO  SUSTAIN  LIMIT  CYCLE  OSCILLATION; 

FLEXIBLE  SURFACE  WITH  PRELOAD  NONLINEARITIES  IN  ROOT  PITCH 
AND  ROLL 


3.3  CORRELATION  WITH  NUMERICAL  RESULTS 


The  asymptotic  solutions  obtained  in  the  preceding  paragraph  were  briefly 
compared  with  the  numerical  simulation  results  of  Reference  1.  The  objective 
here  was  to  compare  the  results  of  the  describing-function  (first  order)  solu¬ 
tions  and  the  asymptotic  (second  order)  solutions  with  the  "exact"  numerical 
solutions  in  the  accuracy  critical  regions.  From  the  results  of  Reference  1 
it  was  concluded  that  the  worst  correlation  between  the  describing-function 
results  and  simulation  results  was  observed  for  the  case  of  a  flexible  control 
surface  with  two  preload  nonlinearities.  The  authors  of  Reference  1  attribu¬ 
ted  this  discrepency  to  the  influence  of  higher  harmonics  neglected  in  the 
describing-function  approach.  For  this  reason  this  configuration  was  used  to 
compare  the  three  solution  approaches. 

The  aeroelastic  response  results  for  the  flexible  baseline  control  sur¬ 
face  with  two  preload  nonl inear i ties  is  shown  in  Figure  21.  The  dynamic  pres¬ 
sure  required  to  sustain  the  limit  cycle  oscillation  is  plotted  as  a  function 
of  the  freeplay-to-amplitude  ratio,  S/A.  The  curves  are  for  a  system  with 
uncoupled  pitch  and  roll  frequencies  of  215Hz  and  140  Hz,  respectively.  In 
Figure  21(a),  the  root  roll  motion  is  shown, and  in  Figure  22(b)  the  root  pitch 
motion  is  shown.  The  results  of  these  comparisons  show  a  better  correlation 
between  the  second  order  asymptotic  solution  and  the  numerical  simulations 
than  was  obtained  using  the  first  order  solutions.  This  supports  the  con¬ 
clusions  of  Reference  1  that  the  source  of  the  error  in  the  one  term  des¬ 
cribing-function  approach  is  a  result  of  the  higher  harmonics  not  being  in¬ 
cluded  in  modeling  the  nonlinear  behavior. 
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FIGURE  21  COMPARISON  OF  PREDICTED  AND  SIMULATION  RESULTS  FOR  A  FLEXIBLE 
CONTROL  SURFACE  WITH  A  PRELOAD  NONLINEARITY  IN  PITCH  AND  ROLL 


4.0  CONCLUSION 


The  presence  of  structural  nonlinearities  can  adversely  affect  the 
aeroelastic  response  of  aerodynamic  surfaces.  Of  particular  concern  is  the 
occurance  of  a  sustained  limit  cycle  response  which  may  lead  to  structural 
damage  of  the  aerodynamic  surface,  support  structure  and  equipment  components 
such  as  control  actuators.  The  results  obtained  during  this  study  show  that 
aerodynamic  surfaces  with  structural  nonlinearities  can  become  susceptible  to 
limit  cycle  behavior  at  dynamic  pressures  below  the  linear  flutter  value.  The 
nature  of  the  limit  cycle  response  was  found  to  be  a  constant  amplitude  motion 
with  a  steady  state  frequency  expressed  as  a  function  of  the  magnitude  of  the 
nonlinearity  and  the  amplitude  of  the  oscillation. 

The  objective  of  the  present  study  was  to  develop  an  analysis  procedure 
to  predict  aerosurface  limit  cycle  response  that  retains  the  flexibility  and 
well  defined  procedures  of  the  describing-function  approach,  yet  provides  greater 
accuracy  and  generality  in  modeling  the  nonlinear  system  behavior.  This  was 
accomplished  by  using  an  asymptotic  expansion  technique  to  derive  a  relationship 
between  the  parameters  characterizing  the  structural  nonlinearity  and  the  amplitude 
and  frequency  of  the  limit  cycle  response. 

The  overall  conclusion  from  this  investigation  is  that  the  use  of  the 
asymptotic  expansion  method  results  in  an  accurate  prediction  of  the  aerodynamic 
surface  limit  cycle  response.  The  procedure  described  in  this  study  was  used  to 
investigate  the  interrelationship  between  the  magnitude  of  the  nonlinearities, 
flight  conditions  and  the  nature  of  the  resulting  limit  cycle  response.  First  and 
second  order  asymptotic  solutions  were  developed  for  a  freeplay  type  nonlinearity 
with  and  without  a  linear  preload  in  the  root  support  structure.  The  asymptotic 
solutions,  which  contain  the  contribution  of  higher  harmonics,  were  compared  with 
the  describing-function  approach  which  contains  only  the  first  harmonic.  It  was 
shown  that  the  first  order  asymptotic  solutions  and  the  describing-function 
approach  were  identical  and  the  higher  order  asymptotic  solutions  could  be  used  to 
account  for  the  contribution  of  higher  harmonics  in  the  limit  cycle  response. 


The  significance  of  the  higher  harmonics  in  predicting  the  nature  of  the 
limit  cycle  response  was  found  to  depend  on  the  particular  nonlinearity  studied. 
During  the  development  of  the  second  order  asymptotic  solutions  for  the  freeplay 
nonlinearity  it  was  concluded  that  the  higher  harmonics  had  small  effect  on  the 
results.  Even  for  amplitudes  near  the  freeplay  magnitudes,  where  the  nonlinear 
effect  is  strongest,  the  higher  harmonics  were  found  to  have  little  influence. 
These  results  are  reasonable  when  one  considers  the  exceptionally  good 
correlation  demonstrated  between  the  first  order,  single  term  solution  and  the 
exact  numerical  solution.  It  is  concluded  that  little  is  gained  in  the  prediction 
of  limit  cycle  response  by  including  the  higher  order  solutions  for  the  freeplay 
nonlinearity. 

In  the  case  of  the  preload  nonlinearity  the  influence  of  higher  harmonics 
was  found  to  have  a  definite  impact  on  the  predicted  limit  cycle  response.  When 
the  first  and  second  order  solutions  were  compared  for  this  nonlinearity,  the 
latter  predicted  a  considerably  greater  reduction  in  the  effective  stiffness  at 
amplitudes  near  the  freeplay  region  of  the  load-displacement  curve.  This  behavior 
was  observed  to  be  consistent  with  the  numerical  simulation  results  presented  in 
Reference  1.  As  the  amplitude  of  motion  became  larger  than  the  preload-plus- 
freeplay  magnitude,  the  first  and  second  order  asymptotic  approximations  of  the 
effective  stiffness  converged.  As  amplitudes  increased  further  both  solutions 
approach  the  linear  solution  results.  In  addition  to  the  influence  on  the 
effective  stiffness  values  the  higher  harmonics  were  also  seen  to  contribute  to 
the  form  of  the  expressions  to  determine  the  limit  cycle  time  history  response. 
Variations  of  amplitude  in  the  computed  waveform  for  the  limit  cycle  response 
could  be  attributed  directly  to  these  higher  order  harmonics.  It  was  concluded 
from  these  results  that  for  the  preload  type  nonlinearity  the  second  order 
asymptotic  solutions  should  be  used  to  accurately  predict  the  nature  of  the  limit 
cycle  response. 

This  study  has  shown  the  applicability  of  the  asymptotic  expansion  approach 
to  account  for  the  Influence  of  structural  nonlinearities  in  the  limit  cycle 
response  analysis  of  aerodynamic  surfaces.  The  method  developed  employs  the 
asymptotic  solutions  to  determine  the  effective  system  parameters  governing  the 
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nonlinear  response.  The  ability  to  include  the  influence  of  higher  harmonics  in 
the  nonlinear  load-displacement  relationship  and  obtain  solution  accuracies  to  any 
desired  order  have  been  demonstrated.  The  methods  were  applied  to  two  simple 
nonlinear  systems  and  the  potential  influence  of  higher  order  solutions  on  the 
limit  cycle  response  was  investigated.  The  applicability  of  the  asymptotic  methods, 
however,  is  not  restricted  to  these  nonlinearities  or  undamped  systems. 
Investigations  of  the  asymptotic  solutions  for  systems  with  other  types  of  non- 
linearities,  including  nonlinear  damping  dependent  on  the  displacement  and/or  its 
derivatives,  are  feasible  using  this  approach. 
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n  General  indicies 

0  Associated  with  root  pitch  degree  of  freedom 

9  Associated  with  root  roll  degree  of  freedom 
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APPENDIX  A 

Flutter  Analysis  of  Baseline  Aerodynamic  Surface 

The  aeroelastic  analysis  performed  in  this  study  utilized  the  flutter 
results  for  the  effective  or  linearized  representation  of  the  baseline  aero¬ 
dynamic  surface.  The  physical  description  of  the  baseline  surface  and  the 
flutter  results  for  both  rigid  and  flexible  cases  are  discussed  below. 

Properties  of  the  Harpoon  missile  control  surface  were  used  to  define  the 
baseline  surface  configuration  which  was  used  throughout  the  study.  The  ge¬ 
ometric  configuration  of  the  control  surface  is  shown  in  Figure  A-l .  The 
structural  nonlinearities  that  were  investigated  are  associated  with  the  root 
support.  Presented  in  Figure  A-2  are  the  inertia  properties  of  the  control 
surface.  The  first  two  rows  and  the  columns  of  the  inertia  matrix  are  asso¬ 
ciated  with  rigid  root  roll  and  pitch  motions  while  the  last  two  diagonal 
elements  are  the  generalized  masses  of  the  control  surface  modes.  The 
off-diagonal  terms,  the  PF  quantities,  represent  the  inertia  coupling  between 
rigid  and  flexible  motions.  The  mode  shapes  associated  with  the  first  two 
control  surface  cantilever  modes  are  given  in  Figure  A-3.  These  modal  data 
were  used  when  investigating  a  flexible  control  surface  configuration. 

Representative  flutter  results  for  the  baseline  control  surface  take  the 
form  illustrated  in  Figures  A-4  and  A-5.  The  results  given  in  Figure  A-4  are 
for  a  rigid  fin  while  those  for  the  flexible  fin  are  given  in  Figure  A-5. 
These  two  figures  show  the  variation  in  the  flutter  critical  dynamic  pressure 
as  a  function  of  the  effective  root  roll  frequency  for  various  values  of  the 
effective  root  roll  frequency. 


*9  *99 

PF9. 

PF92 

P%1 

PF62 

ml 

0 

Symmetric 

m2 

(a)  Form  of  Inertia  matrix 


l0  •  0.1667  lb-see2-ln. 

I.  -  0.071  1b-sec2-ln. 

9 

1b-sec2-ln. 

m^  •  6.220  x  10  4  1b-sec2/In. 
n»2  •  6.295  x  10  ^  1b-sec2/tn. 


PFel  -  -7.227  x  lO-3  lb-sec2 
FFW  -  -1.016  x  10_3  lb-sec2 

P%2  "  “3’212  x  ,0"3  ,b_#*c2 

PF^2  -  2.136  x  10’3  1 b-$oc2 


(b)  Spool  fie  Inert  to  terms 
Figure  A-2  Control  surface  Inert  Is  properties 


FIGURE  A-4  EFFECTIVE  SYSTEM  RIGID  CONTROL  SURFACE  FLUTTER  RESULTS 


ltd  -  imn td  oimtv Aq 


tffectlve  rout  ret  I  frequency  Of  -  Hz 

(b)  AMPLITUDE  RATIO 


FIGURE  A- 5  EFFECTIVE  SYSTEM  FLEXIBLE  CONTROL  SURFACE  FLUTTER  RESULTS 
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APPENDIX  B 

DETAILED  DEVELOPMENT  OF  THE  ASYMPTOTIC  SOLUTIONS 


In  the  discussions  of  the  asymptotic  solutions  presented  in  Section  2.0, 
the  first  and  second  order  approximations  to  the  limit  cycle  response  were 
presented.  The  details  of  the  computational  steps  are  given  here  for  the 
freeplay  and  preload  nonlinearities. 


As  discussed  in  Section  2.0  the  nonlinear  equations  of  motion  for  the 
system  considered  in  this  study  can  be  written  as 


d2x 

7? 


'X  *  e 


f ( x)  +  Q(x) 


(Bl) 


In  the  asymptotic  expansion  methods,  the  solution  of  Equation  (Bl)  is  assumed 
to  be  of  the  form 


N 

X  *  A  COS  <||  +  £ 
n*l 


en  Un  (A,*)  +  0  (eN+1) 


(B2) 


The  amplitude  A,  and  phasing  parameter  iff  are  determined  from 


M  .  y 

rft 


N 


dt 


e"a 


n«l 


+  0  (eN+1) 


di.  w  +  E  en  e_  *  0  UN+1) 

0  n*l  n 


"3T 


(B3) 

(B4) 


The  term  N  indicates  the  order  of  the  asymptotic  approximation.  For  this 
study  solutions  up  to  the  second  order  (N=2)  were  determined.  The /3n  and  Bn 
terms  of  Equations  (B3)  and  (Bl)  are  functions  that  satisfy  Equation  (Bl)  and 
are  expressed  in  terms  of  coefficients  obtained  from  the  Fourier  Series  ex¬ 
pansion  of  the  nonlinear  function  f(X)  in  Equation  (Bl). 


♦In  Appendix  B,  the  number  in  the  second  bracket  identifies  identical 
equations  from  Section  2.0.  As  equations  (B-l),  (B-2),  (B-3)  and  (B-4)  were 
discussed  in  detail  in  Section  2.0  they  will  just  be  presented  here.  Refer  to 
Section  2.0  for  a  complete  description  of  these  equations. 
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The  first  step  in  the  solution  procedure  is  to  determine  functions  o^, 
a2’  61 *  e2*  and  ^hese  functions  can  be  derived  by  expanding  the  nonlinear 
function  f  (ip )  in  a  Fourier  Series,  substituting  Equations  82  through  B4  Into 
Equation  81  and  equating  coefficients  of  sin  $  and  cos  if».  These  functions  are 
given  by  the  following  relations. 


*1  2ttu  A 


f  < 

o 


f  (A  cos  i|i)  sin  ipdip 


f  (A  cos  <t»)  cos  <(>dtj> 


n  .  m  g  cos  n*  +  h  sin  n ip 

U,  ■  -4  -  -U  £  - T— 2 - 

'  V  “o2  n-2  n-1 


where 


1  C2lt 

g r\  *  Tri  J  f  (A  cos  *)  cos  o'l'difr  (E 

i  rZlT 

hn  *  lit  J  *  (A  cos  ^  sin  n,|'d',;  (E 

The  terms  «],  3],  gn  and  hn  are  the  coefficients  of  the  Fourier  Series  ex¬ 
pansion  of  the  nonlinear  force  f ( X )  acting  on  the  system.  The  second  order 
coefficients  <*2  and  3 2  can  be  expressed  in  terms  of  Equations  (B5)  through 
(B7)  and  the  derivative  of  f(X)  as 


“2  "i:  |2ai  si  **!  is)n*d» 

1  i  2  ®1  ^*1  1  (  f2* 

b2  "^j8i  *r  ±cos 

I  ■  ^  f'%  (a  cos  1|>)  ♦  cos  *  -  Ab1  sin  *  +  U  f  J 


(BIO) 


(Bll) 


(B12) 
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By  computing  the  terms  in  Equations  (B5)  through  (B12),  substituting  into 
Equations  (B3)  and  (B4)  and  integrating,  the  amplitude  A  and  phasing 
parameter,  tfi  approximated  to  the  second  order,  are  obtained. 


Equations  (B5)  through  (812)  are  for  the  general  case.  For  the  two  non- 
linearities  considered  in  this  study,  these  equations  can  be  reduced  to  a  much 
simpler  form.  First,  the  form  of  f ( X)  for  the  freeplay  and  preload  nonlin¬ 
earities  is  such  that  all  <*n  and  hn  terms  in  Equations  (B5),  (BIO),  and  (B9) 
are  identically  zero.  For  this  case.  Equations  (B3)  and  (B4)  yield 


A  *  const 


*  * 

or 


ut 


CJ 


(u»o  +  cB1 

f 

u)  +  eBi 

l  o  1 


+  e2  B2jt  +  const. 

*  '2  h) 


(B13) 

(B14) 

(B15) 


where  55  is  the  effective  frequency  approximated  to  the  second  order.  Equa¬ 
tions  (B7)  and  (Bl)  then  become 


U 


1 


0» 

L 

n=2 


gn  cos  * 
n2-l 


(B16) 


(B17) 


The  computational  effort  of  this  study  was  centered  on  determining  Ut ,  3]  and 
029iven  by  Equations  (B6),  (B16)  and  ( Bl 7 ) ,  respectively. 


The  form  of  the  effective  stiffness  to  the  second  order  given  by  Equa¬ 
tion  (4)  is 


7  •  K 


^1  +  c  B-j  +  e2  B2  +  •  •  •  c 


( Bl  8) 
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This  expression  can  be  obtained  from  Equation  (B1 5)  by  squaring  both  sides  and 
retaining  term  up  to  <  2.  This  yields. 


(B19) 


Substituting  Equations  (B6)  and  (B17)  for  and  02  into  Equation  (B19)  and 
simplifying  yields 


T3^  ■  «o  £l  +  c  Bi  +  e2  Bjj 

The  relationship  between  B]  and  B2  and  £]  is  defined  as 


(B20) 


(B21 ) 
(B22) 


where  <J>  is  defined  by  equation  B-12. 
Note  that  here  it  is  implied  that 


7  /*2ir 

B2  *  ■  izn  J  £ cos 

o  0 

. .  _2  2 
lim  w  *  w 

0 

c^o 


(B23) 

(B24) 


Using  the  relationship  between  the  effective  stiffness  and  natural  frequency. 
Equation  (B18)  is  obtained  directly  from  (B20). 


Summarizing  the  computational  procedure,  the  first  step  is  to  compute  0 
for  the  linear  system.  Then  the  term  fi-\  is  obtained  from  Equation  (B6).  Next 
coefficients  Bi  and  B2  can  be  determined  using  Equations  (B12),  (B21)  and 
(B22).  At  this  point  the  second  order  correction  to  the  effective  stiffness 
is  obtained  from  Equation  (B18).  The  computation  of  the  time  domain  solution 
then  requires  solution  of  Equations  (B8),  (B14)  and  (B16).  The  form  of  these 
relationships  are  presented  below  for  the  two  nonlinearities  considered  in 
this  study. 
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FREEPLAY  NONLLINEARITY 


In  order  to  compute  the  coefficients  of  the  asymptotic  approximations, 
the  initial  step  is  to  expand  the  load  relationships,  l(t)  in  a  Fourier 
Series.  This  load  relationship  is  shown  in  Figure  4  of  the  text  and  is 
expressed  as 


L(t) 


Io  Ifor  0  <  t  < 

K[A  sin  t  -  S](for  t,  <  t  <  2ir 


(B25) 


This  load  is  in  the  form  of  the  function  F(Acos^)  in  Equation  (B21)  so  the 
first  Fourier  coefficient  corresponds  to  the  Bi  term  in  Equation  (B20).  Since 
the  load  function  l{t)  is  a  single-valued  odd  function  its  Fourier  Series 
representation  is  of  the  form 


b  cos  n  t 
n 


(B26) 


/2tt 

i  (t)  cos  0  t 

o 


(B27) 


As  mentioned  above  the  b-j  term  of  Equation  (B27)  and  B]  in  Equation  (B20)  are 
related  by  the  expression 


'  ■(*=? b’  - ') 


(B28) 


The  form  of  the  b]  term  for  the  load  of  Equation  (B25)  is, 
A u  2  /  \ 

bi  ’  ir~  v  * 2  h  m  sin  2  V 


(B29) 


Computation  of  the  B2  term  of  Equation  (B23)  requires  the  coefficients  of  the 
higher  harmonics  of  Equation  (B26).  A  general  expression  for  these  higher 
harmonics  is  determined  as 


b  *  0  (n  ■  2,4,...) 
n 


b  •  !jL-  Fp  -  -  sin  n  t,1  (n  *  1,3,5...) 

n  it  L  N  n  J 


(B30) 
(B31 ) 


(B32) 


r 


n 


sin  (n-1)^  sin  (n+l)t.,’ 
2  (n-rr  +  2  (n+1 )  . 


By  successive  substitution  of  Equation  (B30)  into  Equations  (B16),  (B12) 
and  (823)  the  82  coefficient  of  Equation  (B18)  is  obtained.  The  form  for  the 
freeplay  nonlinearity  is  given  in  Equation  (B33). 


1  E 


n»3.5 


t^  sin  nt 


1 


(B33) 


PRELOAO  NONLINEARITY 


The  procedure  for  deriving  the  coefficients  of  the  asymptotic  approxi¬ 
mation  of  Equation  (B18)  for  the  preload  nonlinearity  follows  the  same  proce¬ 
dure  used  for  the  freeplay  nonlinearity.  The  load -displacement  relation  for 
the  preload  nonlinearity  is  shown  in  Figure  7  and  is  defined  as 


I(Aq  ♦  A1  cos  t  -  2S)K 
PK 

(Aq  ♦  A1  cos  t)K 
PK 

(Aq  +  A1  cos  t  -  2S)K 


for  0  <  t  <  tj 

for  t,  <  t  <  t2  (B34) 

for  tj  <  t  <  2t  -  t2 
for  2*  -  t2  <  t  <  2*  -  t| 
for  2»  -  t^  <  t  <  2ir 


t1  «  cos 


P  *  ZS  -  A„ 

( - *> 


(B35) 


1  P  '  A« 
COS  1  (  -A-—) 
A1 


The  Fourier  Series  expansion  of  the  load  relationship  is  defined  as 

b  “ 

L  (t)  ■  y-  +  £  bn  cos  nt  +  «n  sin  nt 
n«l 


(B36) 


(B37) 


bQ  *  y  L  (t)  dt  (B38) 

bn  is  given  by  Equation  (B27)  and  all  <*n  coefficients  are  zero  for  the  func¬ 
tion  given  by  Equation  (B34). 


The  bn  coefficients  of  Equation  (B35)  for  the  preload  nonlinearity  where 
P  <  A  <  P+2S  are 


bo  ■  -V1  [2  Ao  (*-‘1 )  ♦  2  Pt,  -2A,  »<»  t,] 

51  ’  -V1  [’-»!  *  *7  (p-A0)  s'"  *,  -  7  sfn  2  t,j 

•  2A,  \l  (P-A  )  1 


(B39) 

(B40) 

(B41) 


t]  is  defined  in  Equation  (B35)  and  rn  is  defined  in  Equation  (B32).  For 
amplitude  of  motion  A>P+2S 


u  2A 

bo  =  -2“ 
0  IT 


ui  2a 

‘l  '-f- 


*7  (t2+1f*tl)+  (s,n  *2  •  si"  tl)  *  *7  (tT*z)-  A*  *< 


(B42) 


’rt2'tl  '  I7  (P*2S-Ao)  s,n  *2 
*  17  (P-Ao)  sin  ll  *  I  (sin  2  1 


2  -  sin  2  t1 


(B43) 


„2a 


(Ao-p)(sin  nt^-sin  nt^  -2S  ^sin  nt^  +  on 
t,  and  t2  are  defined  by  Equations  (B35)  and  (B36)  and  <rn  is  defined  as 


(B44) 


A 


sin  (n-l)t2  -  sin  (n-1 ) t-j  sin  (n+l)t2  -  sin  (n+l)^ 


"nrrrr 


tt^w 


The  detailed  derivation  of  the  asymptotic  solutions  for  both  the  freeplay 
and  preload  nonlinearities  have  been  presented  here.  These  derivations  were 
included  to  supplement  the  discussions  of  section  2.0  and  should  be  referenced 
when  more  detail  is  required. 
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